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Three-dimensional topological insulators are characterized by the presence of protected gapless 
spin helical surface states. In realistic samples these surface states are extended from one surface 
to another, covering the entire sample. Generally, on a curved surface of a topological insulator an 
electron in a surface state acquires a spin Berry phase as an expression of the constraint that the 
effective surface spin must follow the tangential surface of real space geometry. Such a Berry phase 
adds up to 7r when the electron encircles, e.g., once around a cylinder. Realistic topological insulators 
compounds are also often layered, i.e., are anisotropic. We demonstrate explicitly the existence of 
such a 7r Berry phase in the presence and absence (due to crystal anisotropy) of cylindrical symmetry, 
that is, regardless of fulfilling the spin-to-surface locking condition. The robustness of the spin Berry 
phase 7r against cylindrical symmetry breaking is confirmed numerically using a tight-binding model 
implementation of a topological insulator nanowire penetrated by a 7r-flux tube. 



I. INTRODUCTION 

Consider the situation where the lower half of a three 
dimensional space is occupied by a topological insulator 
with Z 2 index v = 1 and the rest is a vacuum (y = 0). 
Then as the most characteristic feature of the topological 
insulating state, a metallic surface state appears at its in- 
terface with the vacuumP^ An interesting variant of this 
scenario explored by a recent Aharonov-Bohm measure- 
ment on Bi2Se3 nanowire^ and subsequent theoretical 
analyses^^l is that the surface metallic state is not only 
protected by time-reversal symmetry, but it shows an- 
other characteristic feature when the surface is deformed, 
say, into a cylinder (see FIG. 1) — the manifestation of 
the spin Berry phase. Depending on how the surface is 
deformed into a cylinder, i.e., whether the topological 
material fills either the inside or the outside of the cylin- 
der, the system can be regarded as either a nanowire 
or a linear aperture penetrating an otherwise surfaceless 
topological insulator. In a recent worlf^l we have cho- 
sen the latter as a starting point for studying the nature 
of topol ogicall y protected helical modes along a disloca- 
tion line! 11 * 12 ! Effects of a finite size (a finite radius of 
the cylinder) combined with the presence of a nontrivial 
spin Berry phase was shown to play an essential role in 
protecting the ID helical modes JSnZUQI 

The appearance of spin Berry phase n is a character- 
istic feature of topological insulator surface stated dis- 
tinguishing it from, e.g., a carbon nanotubeJSl another 
2D gapless Dirac system (i.e., graphene) rolled up into a 
cylinder. These two Dirac systems both involve an effec- 
tive spin degree of freedom appearing in the low-energy 
effective Hamiltonian. The physical origins of these spin 
degrees of freedom are different in the two cases; in the 
case of graphene (or carbon nanotube) it is the sub-lattice 
structure of hexagonal lattice, whereas in the present ex- 
ample it is essentially a genuine electron spin. These two 
effective spin degrees of freedom, despite their very differ- 
ent nature, play a similar role in determining the trans- 
port characteristics of the surface states on a flat surface. 




FIG. 1. Cylindrical surface of a topological insulator nanowire 
and a typical spin configuration of a surface electronic state 
(left). Top view (a cross section is shown) of the spin configu- 
ration (right) in the (a) presence and (b) absence of cylindrical 
symmetry. In case (b) the spin follows the tangential plane 
of an auxiliary elliptic surface. 



This is, however, no longer the case when the surface is 
curved. The sub-lattice pseudo-spin of graphene is in- 
sensitive to warping of the 2D plane. The effective spin 
on the surface of a topological insulator is, on the con- 
trary, constrained to lie in-plane to the surface.^ This 
constraint is the origin of the spin Berry phase it charac- 
teristic to the topological insulator surface statePtH One 
of the purposes of the present work is to demonstrate 
through explicit examples how the information encoded 
within the bulk Hamiltonian manifests itself in the sur- 
face effective Hamiltonian in the form of a nontrivial spin 
Berry phase. In the course of deriving the spin Berry 
phase, we also establish an unambiguous correspondence 
between the effective spin degree of freedom appearing in 
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the 2D surface Dirac Hamiltonian and the original real 
spin embedded into the 3D bulk effective Hamiltonian. 

A second motivation of this work is to explore the con- 
sequences of the anisotropy of topological insulators on 
the surface states, especially on the spin Berry phase. 
A rectangular nanowire made of such asymmetric com- 
pounds have surfaces of different symmetries. For ex- 
ample, in the Bi2Se3 nanowire studied in RefP the sur- 
faces orthogonal to the c-axis is rotationally symmetric, 
whereas surfaces parallel to the c-axis is not. Correspond- 
ingly, the surface Dirac cones are symmetric in the for- 
mer, whereas distorted in the latter. The slope of the 
energy dispersion (Fermi velocity) also differs. What 
happens to the electron when it goes through (or gets 
reflected? by) junctions between two surfaces of different 
character? What is the fate of the Berry phase tt when 
an electron goes around the nanowire passing by several 
of such junctions? Such effects of anisotropy and multi- 
ple surface geometry will be also important in a similar 
Josephson-Majorana geometry involving metallic surface 
and Majorana bound statesP^l 

This paper is organized as follows: we first demon- 
strate (Sec. II) by a simple analytic calculation that in 
the presence of anisotropy the surface effective spin is 
no longer strictly locked in-plane to the tangential plane 
on a curved surface, but it has generally a finite out-of- 
plane componentP21 We also show that only the global 
Berry phase tt which an electron acquires when it winds 
around a cylinder is robust against such asymmetry of 
the crystal. This point is further confirmed numerically 
by implementing a topological insulator nanowire pen- 
etrated by a 7r-flux tube as a tight-binding model on a 
square lattice (Sec. III). 



II. DERIVATION OF THE SPIN BERRY PHASE 

Let us first derive the spin Berry phase tt directly from 
a bulk 3D effective Hamiltonian. In contrast to RefsP^, 
here, we choose to go back to the 3D bulk effective Hamil- 
tonian, and derive the spin Berry phase directly from the 
3D Hamiltonian. 



Cylindrical nanowire in parallel with the crystal 
growth axis 



Let us start with the case in which the system has a 
cylindrical symmetry, i.e., a bulk insulating state is con- 
fined inside a cylinder of a radius R: \J x 1 + y 2 < R, 
directed along a crystal growth axis (c-axis) perpen- 
dicular to the stacking layers. To describe the bulk 
insulating state, which is in contact with the vacuum 
outside the cylinder, we take the following effective 



HamiltonianpSHin] 



H3D = 



M Bk z Ak,_ 

Bk z -M Ak- 

Ak+ M -Bk z 

Ak, -Bk 7 -M 



(1) 



This Hamiltonian describes a 3D Dirac system with a 
mass parameter M = Al + M 2 (k x + k 2 + k 2 ) . The nature 
of the insulating state, i.e., whether it is Z2 trivial or not, 
is determined by the sign of M /M 2 . When M Q /M 2 < 0, 
the insulating state is Z2 nontrivial: v = 1 (the gap is 
inverted) exhibiting a single gapless surface Dirac cone, 
whereas when M$/M 2 < 0, Z 2 index is v = and the 
gap is normal. We regard here A and B to be constant; 
A = Aq, B = Bo in the parametrization of Rcf. 19 . Note 
that the same is assumed in the derivation of flat surface 
states in RefpS. The anisotropy of the crystal is reflected 
in the asymmetry between A and B; according to Table 
IV of Ref.™, Bq is generally smaller then Ao, and in the 
case of Bi2Te3 smaller by one order of magnitude. The 
structure of the Hamiltonian ([!]) may become clearer in 
the follwoing symbolic form: 

H = Mt z + Bk z T x a z + Ar x (k x a x + k y a y ), (2) 

where r = (t x ,t v ,t z ) and er = (a x , a y ,a z ) are two sets of 
Pauli matrices representing, respectively an orbital and 
a spin degrees of freedom. 

To identify the gapless surface states in the cylin drical 
geometry, we first decompose, in the spirit of Ref) 19 * 20 ^, 
the bulk 3D effective Hamiltonian ([lj into two parts: 



H = H ± (k r ) + H {] (k^,k z ), 



(3) 



where H± = H\k^=k z =o- k r and fc^ are components of 
the crystal momentum conjugate to the cylindrical coor- 
dinates: 



V x2 + 



y 



y 

arctan — . 

x 



(4) 



H± is a Hamiltonian at the T-point, and reads explicitly, 



H±_ = M\t z + Ak r T x (a x cos </> + a y sin </>) (5) 
Mi. Ae-^k, 

-Ml Ae-^K 
Ae^K Mi 
Ae^kr — Mj_ 



(6) 



For a later use, note also that the remaining H\\ becomes 



My 
Bk z 


iAe^kd 



Bk z 
-My 
iAe^k^ 






-iAe~ 

-Bk 



*k, 



-iAe- % 


-Bk 
-Mi 



*k d 



(7) 

In Eqs. ^ and Q we have decomposed the mass term 
into Mi ~ M +M 2 k 2 and M\\ = M 2 {k\+k 2 z ) . Of course, 
the Laplacian in the cylindrical coordinates has another 
contribution, (\ / r)d / dr . Here, we neglect this first-order 
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derivative term, keeping the term d 2 /dr 2 , which is a pos- 
teori justified, since the penetration depth A of the sur- 
face state is much smaller than the radius of the cylinder; 
R > A. 

We then consider a solution of the eigenvalue equation, 



H±\il>) =E ± \i/>) 



of the form!™] 



IV-) 



(8) 



(9) 



i.e., k r = —iX (we keep only A > 0). For a given E±, one 
finds four independent solutions of this form, \ip) = \ipj) 
(j = 1, 2, 3, 4). Then one composes a linear combination 
of these four solutions, 



3=1 

for satisfying the boundary condition: 



(10) 



lr=R 



(11) 



The boundary condition (111 gives a restriction to the 
spinor part of the wave function \ip), allowing for explic- 
itly writing down its two independent bases. 

To find the eigenspinors, which compose the spinor 
part of one needs to diagonalize H±. In view of 

its specific form written symbolically as in Eq. |5]), one 
may diagonalize its real spin part first. Namely, one can 
partially diagonalize H±_ as, 



H ± \r±) = (M x t z ± Ak r r x ) \r±), 



(12) 



in terms of the following (real) spin eigenstates pointed 
in the radial direction f = (cos <f>, sin 4>); 



\r+) 



1 

7s 



o»0 



1 

7^ 



Z(f> 



(13) 



Here, we have chosen these eigenspinors single-valued. It 
is possible, of course, to take them double valued, but the 
two choices turn out to be completely equivalent. The 
advantage of the single-valued choice is that the Berry 
phase becomes explicit in the surface effective Hamilto- 



nian; see Eq. (23 1. Whether one chooses one set of |r±) 
or the other, they compose a set of bases diagonalizing 
real spin states pointed in the direction of r. 

The remaining orbital (t-) part can be also diagonal- 
ized as, 



H ± \X±)\r±) =E ± \X±)\r±), 

1 



|A±) = 



±i(E ± -M ± )/(XA) 



(14) 



In the second line we took k r = —iX explicitly into ac- 
count. For a given energy E± satisfying E\ = Mq — 



M%\ , there are four possible solutions for the penetra- 
tion depth, A = iA^', of which we keep only the two 
positive solutions, A = A^. For each of A = A^ we 
have two independent base spinors |A±); we have in to- 
tal four independent solutions, constituting the general 
solution ( |10[ ). To be explicit, the general solution (10) 
reads explicitly, 

\1>) = [ci|A< + >+)|r+)+ca|A< + >->|r-)] e A(+)(r ^ } 
+ [c 3 |A ( - ) +)|r+) + c 4 \X^ -)\r-)] e x ^ ^ R \l5) 



Notice that in Eq. 
of A; E± = E±(\W 



( 14 ) E± and M± are also functions 



E\ ±} , Mi 



M 



(±) 



The last step is to impose the boundary condition ( j'l 1 [ ) 
to Eq. §1^. 

Since the two base spinors \r+) and \r+) subtend or- 
thogonal real spin subspaces, the boundary condition 
(111 requires that 



ci|A«- 
c 2 |A«- 



•c 3 |Af-)- 
c 4 |A(-)- 



(16) 



independently hold. This means that the spinor part of 
Eq. (15) can be expressed solely in terms of the eigen 
spinors, say, with A = A^ + -*, i.e. 
of|AW+) and|A(+>- 



as a linear combination 



det 



iX^A 



E { + ] - M^' E]_> -M ± 



The first line of Eq. ( |16[ ) implies, 

(17) 



i\(-)A 

(-) _ A#(-) 



0. 



One can verify that this holds true only when the two 
conditions: (i) MqM^ < (the system is in the v = 
1 phase) and (ii) E± = arc simultaneously satisfied. 
Substitute E± — into Eq. ( 113 ), and notice that Mj_ — 
—XA is the only choice consistent with the requirement 
that A > 0, if one defines the parameters such that A > 
and M 2 > 0. 

Taking all these into account one can express the so- 
lution of the boundary problem as, 



Cl 



|A( +) +)|r+> + C2 |A(+)-)|r-)l p(r) (18) 
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ie l<t > 




ie lct > 





p{r) 



where 



p{r) 



Ci|r+))+ca|r-» 



V /A+A_(A+ + X-)/nR 



(19) 



|A+-A_| 
(A± <C R assumed) and, 



e \+(r-R) _ e A-(r-R) 



Ah 



A ± y/A 2 + 4M M 2 
4Mo ' 



(20) 



(21) 
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The four-component eigenspinors |r±)) introduced in 
Eq. ( 19 1 describe electronic states localized in the vicin- 



ity of the surface. The effective 2D surface Hamiltonian 
H2D is obtained by calculating the matrix elements of 
Hn in terms of these |r±)), i.e., 



(H 2D ) ±± = ((r±\H l] \r±)). 
By an explicit calculation, Hid is found to be 
H2D = 



(22) 



(23) 







-iBk 7 



a 



a ( 

r\ 1 a<p 





where we have used, k^ ~ —i(l/R)d/d(f>, since R ^> A, 
i.e., only r ~ R is relevant. Two factors 1/2 which have 
appeared in the off-diagonal elements of Eq. ( 23 1 are the 



spin Berry phase terms, which lead to a 7r-phase shift 
when an electron goes around the cylinder. This Berry 
phase term can be eliminated from the eigenvalue equa- 
tion for H2D 



H 



2D 



Cl 


= En 


Cl 


C2 




C2 



by introducing a singular gauge transformation, 



c = 


Cl 


= e -i^/2 


Xl 




c 2 _ 




X2 



(24) 



(25) 



In the transformed x-basis, the surface Hamiltonian takes 
a simple Dirac form without the Berry phase, 



H. 



(x) 



Bo y k z , 



L 2D ~ Ao~ x k, 
whereas the corresponding eigenspinors, 



X2 



(26) 



(27) 



become double- valued. 

It is suggestive to express x explicitly in terms of a set 
of polar coordinates defined in the effective surface spin 
space. By introducnig the parameters as, 



cos 77 — 



Ak, 



A 2 kl + Bk 2 



sin 77 = 



BK 



A 2 k 2 +Bk 2 



(28) 



the eigenstates of Eq. ( 26 1, corresponding to the eigenen- 
ergies, 



Eu=±^A*kl + Bk 2 z , 
can be written, respectively, as 



1 

71 



e -«7/2 



(29) 



(30) 




sy 



T 

sx 



FIG. 2. Spin coordinates on the cylindrical surface. 

Comparing this with the textbook formula of an SU(2) 
spinor, 



\n±) 



e -i V /2 cos (Q/2) 

ie^/ 2 sin(6»/2) 



(31) 



pointed in the direction of a unit vector n specified by 
a polar angle 9 and an azimuthal angle 77 of the SU(2) 
spin space: (s x ,s y ,s z ), one can verify that the surface 
effective spin is locked in the (s x , Sj,)-plane, i.e., 9 = tt/2. 

In the derivation of Eq. (23), we have chosen the spin 
quantization axis in the direction of f; see Eq. ( |13[ ). 
Then, on which plane is the surface (ci, C2)-spin actually 
locked? In accordance with Eq. ( p~3| ) , the spin-space 
coordinates should be redefined as 



Sz 



X — > —z, 

0, 



V 

z 



(32) 



on the cylindrical surface (see FIG. 2), since k^ (k z ) play 
the role of k y (—k x ), where 4> = (— sin 0, cos 0). Taking 



this into account, one can interpret Eq. (26), neglecting 
the anisotropy (B 



H. 



A), as 

(x) ~ A{ax k) 



2D 



(33) 



In this regard, one can view Eq. ( |19[ ) with coefficients 
X\ and X2 given in Eq. p0]) as an SU(2) spin state, 



(34) 



\r±) = Xi\r+) + X2\r- 



decorated by an accompanying orbital degree of freedom. 
Eq. (34) represents a spin state which is locked in the 



plane perpendicular to f; unit normal vector of the cylin- 
drical surface. This is a clear indication that the effective 
spin on the surface of a topological insulator cylinder is 
a real spin which is locked at each point of the cylinder 
in-plane to its tangential surface [FIG. 1, panel (a)], un- 
like the sub-lattice pseudo-spin on the cylindrical surface 
of a carbon nanotubeP^ 
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One step backward, notice that this peculiar property 
of the spin Berry phase 7r manifest in Eqs. (23) and 



( 25 ) is here not derived from the property of the effective 

( 26 1 on the cylindrical coor- 



Eq. 



Dirac Hamiltonian, i.e. 
dinates. It was encoded in the dependence of eigenspinors 
( 18 ) on the spatial angle cf>, and the explicit form of H\\ 
as given in Eq. @. In RefsPHthe same conclusion was 
drawn by observing the surface effective Hamiltonian on 
a curved surface. Here, the Dirac equation was derived 
simultaneously with the spin Berry phase; the existence 
of the spin Berry phase is indeed encoded in the bulk 3D 
Hamiltonian. 



B. Cylindrical nanowire perpendicular to the 
crystal growth axis 

We have considered so far electronic states on the sur- 
face of a cylinder whose axis of symmetry is pointed along 
the c-axis, i.e., in the direction of crystal anisotropy. 
Since the rotational symmetry around the c-axis is pre- 
sumed, any tangential surface of the cylinder is equiva- 
lent. Under such circumstances, we have seen explicitly 
that the effective spin degree of freedom appearing in the 
Dirac equation (23) is constrained to the curved (cylin- 
drical) surface. 

Here, we consider a less trivial case of broken rotational 
symmetry, i.e., with the cylindrical axis chosen perpen- 
dicular to the c-axis. Of course, nanowires are not cylin- 
drical in real samples, but have several surfaces. It is also 
unlikely that the axis of the wire is perfectly aligned with 
the axis of crystal symmetry (c-axis). In the presence of 
such rotational anisotropy, it is less trivial whether the 
effective spin degree of freedom on the surface is always 
tied to the curved surface. 

To implement an anisotropic cylindrical surface we 
consider here the case of the crystal c-axis pointed in 
the ir-direction, keeping the symmetry axis of the cylin- 
der pointed always in the z-direction; obviously, one can 
equally rotate the cylinder in the y-direction with keep- 
ing the crystal growth axis in the z-direction. The bulk 
effective Hamiltonian as Eq. ^ for such a rotated crystal 
becomes (here, we do rotate the crystal), 



H 3D = Mt z 



5 {Ba x k x 



Ao~yky + Aa z k z ). (35) 



To identify the surface electronic states which span the 
basis for the 2D surface Dirac Hamiltonian, we introduce 
the same cylindrical coordinate as Eq. Q and decom- 
pose the Hamiltonian (35) into perpendicular (H±) and 
parallel (ify) components in parallel with Eqs. (|3]), (|6j) 
and (JTJ). Some parameters need, of course, redefinition 
or exchange. Let us first focus on 

H± = M\r z + k r r x (Ba x cos <p + Aa v sin <ft) . (36) 



To diagonalize the real spin part of the Hamiltonian ( 36 1 , 



B cos 4> = A cos q 
A sin <fi — A sin ( 



where 



A = A(<j>) = JB 2 cos 2 <j> + A 2 sin 2 , 



(37) 



(38) 



In analogy with Eq. (12), one can partially diagonalize 
Eq. Q as 



H\ \f±) = M 



\f±) 



(39) 



but here the eigenspinors represent no longer spin states 
in the f-direction normal to the surface of the cylin- 
der. The new eigenspino rs \ r±) are formally analogous to 
|r±) defined as in Eqs. (131, but pointed in the direction 
specified by <fr introduced above; 



1 

75 



1 



1 

7^ 



i 



(40) 



The remaining procedure is perfectly in parallel with the 
previous case. After imposing the boundary condition on 
the cylindrical surface of the topological insulator, one 
finds as the basis spinors for constructing the surface ef- 
fective Hamiltonian, 



Cl 



|A<+>+)|f+) + C2|A<+>-)|f-)l p (r) (41) 
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p{r) 



c 1 \f+))+c 2 \f-)}, 



(42) 



where p(r) is given as in Eq. (20) with a normalization 
factor expressed in terms of X± given as in Eq. (21 ) but 
with A replaced by A given in Eq. ( 38 1 . 



To find the surface effective Hamiltonian, we calculate 
again the matrix elements of H\\, here, in terms of |r±)) 
associated with |f±) introduced in Eqs. (40). We have 
introduced the parameter (f> in Eqs. (37) and was able to 



construct the basis spinors which successfully span 

the subspace representing the solution of the boundary 
problem. 0-parametrization was useful, since H± shows 
a nice transformation property in terms of <f>. It is, how- 
ever, no longer the case for H\\ , which reads explicitly in 
the present case, 



it is convenient to introduce an auxiliary angle (f>, defined 



#H = .\/ - • Ar,a z k z 

+ t x {—Ba x sin <p + Ao~ y cos < 

Mm Ak z —ic^k^ 

Ak z —-Mil —ic<f,k,p 

ic*^ My -Ak z 

iAks -Ak z -Mm 



(43) 
(44) 



G 



where we have introduced 



= A cos (f> — iB sin ( 



Let us come back to the question, "In which direction 
is the eigenspinor c, or equivalently x, pointed in the spin 
(45) space?" Taking it into account that our surface effective 



is a complex number which is a function of the real 
parameter cf>. Note that does not simplify with the use 
of 4>. 

To derive the surface effective Hamiltonian H2D, one 
calculates again the matrix elements, ((f ± |Uy|f±)) in 
terms of the new r-basis. One can explicitly verify that 
diagonal elements of the surface effective Hamiltonian 
vanish, i.e., 



H- 



2D 



«f + |JT|||f-» 

({f-\m\f+)) 



(46) 



The off-diagonal elements involve a first-order derivative 
with respect to 0, yielding a Berry phase term; e.g., 



({f+\H {] \f-}) = 

-iAk z + I (c^e^ + e-^cifc, 



(47) 



Note that terms coming from the ^-dependence of A±, 
given as in Eq. (21 ), cancel when integrated over r. The 



^-dependence exists not only in the explicit r-dependence 
of p(r) on A±; cf. Eq. (201, but also in the normalization 
factor of p{r). The ^-dependence of X± stems from that 
of A = A(4>). In Eq. (471 the Berry phase term appears 

in the second part when k^ acts on the exponent of e l< ^, 
and is found to be of the form, 



fi(0) 



c^e 1 ^ d<p c<f,e 



i<t> AB 



2R d(j) 2R A 2 



(48) 



Notice that only the real part of this factor influences 
the phase of the wave function and is identified as the 
local spin Berry phase. This phase factor can be elim- 
inated from the eigenvalue equation: B-idC- = -^|| c i by 
empl oyin g a singular gauge transformation analogous to 
Eq. (25), i.e., by c = e~ l<t3 / 2 X- Since </> and 4> have the 



same winding property, 



2tt) 



2tt, 



(49) 



the eigenspinor, Eq. \21\ is precisely double-valued with 
respect to a 27r-rotation of 4>, i.e., the spin Berry phase 

iS juSt 7T. 



The off-diagonal elements of Eq. ( 46 1 are also suscep- 
tible of the imaginary part of £l(<fr) given as in Eq. (48). 
This imaginary part introduces a non-uniform modula- 
tion into the amplitude of the wave function of surface 
electronic states. Such a situation is somewhat similar to 
the case of a WKB wave function describing the system 
with a spatial non-uniformity. The phase of such WKB 
wave function is expressed in terms of a "wave vector" 
k which depends on the space coordinate. Meanwhile its 
amplitude shows also a space dependence. In the present 
formulation, this appears as an imaginary part of Q(4>). 



Hamiltonian is defined on a space spanned by |^±)), we 
redefine the spin coordinates as, 



&x ■ 


x - 


■> — z 


Sy : 


y - 




Sz ■ 


z - 


> r, 



(50) 

where <f) = (— sin0, cos</>). T his means that the eigen- 
state c of H-2D given in Eq. ( 46 1 can be viewed as an 
SU(2) spin state, 



\f±) = bx\f- 



b 2 \r-) 



(51) 



( 51 ) represents a spin 
Note that f 
Thus 



accompanied by orbital spins. Eq. 
state locked in the plane perpendicular to f. 
is not a unit vector normal to the cylinder surface, 
in the case of broken cylindrical symmetry, the surface ef- 
fective spin can have a finite component normal to a tan- 
gential surface of the cylinder [FIG. 1, panel (b)]. An out- 
of-plane component of the surface effective spin appears 
also as a consequence of the hexagonal warpingp2t23 and 
in a system of topological insulator quantum dotP^ 

In this section, we have examined the electronic states 
on cylindrical surfaces of an anisotropic topological insu- 
lator. An explicit one-to-one correspondence between the 
effective spin in the surface Dirac Hamiltonian and the 
real spin inherent to the 3D bulk effective Hamiltonian 
has been established. The existence of spin Berry phase 
7r has been unambiguously shown in this context. In the 
cylindrically symmetric case (Sec. II A), i.e., in the case 
of a cylindrical nanowire parallel to the crystal c-axis, 
we have shown explicitly that the effective surface spin is 
constrained in-plane to the real-space tangential plane of 
the cylinder. In the absence of such cylindrical symmetry 
(Sec. II B: case of a cylindrical nanowire perpendicular 
to the crystal c-axis) the effective surface spin can have 
a finite amplitude in the direction normal to the cylin- 
drical surface. However, when the reference point on the 
cylindrical surface travels once around the cylinder (i.e., 
winds the cylinder once), the resulting spin Berry phase 
is indeed it. This point will be further confirmed in the 
numerical experiments in Sec. III. 



III. ENERGY SPECTRUM IN THE PRESENCE 
OF A tt-FLUX TUBE 

The spin Berry pha se a ppearing in the surface Dirac 
Hamiltonian; cf . Eq. ( 23 ) , is often discussecP^ in the 



context of Aharonov-Bohm experiment on topological in- 
sulator nanowiresF' Let us first recall that the existence 
of a spin Berry phase ir leads, on the surface of a cylin- 
drical nanowire, to the appearance of a finite-size energy 
gap in the spectrum of surface electronic states. The spin 
Berry phase tt modifies the periodic boundary condition 
around the wire to an anti-periodic boundary condition, 
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FIG. 3. Energy spectrum E(k z ) of a topological insulator 
nanowire with a rectangular cross section (N x = 12, N y = 
16). The horizontal axis k z represents the crystal momentum 
along the wire. The crystal anisotropy is introduced in the 
direction perpendicular to the wire (in the x-direction, B/A = 
0.7). M /M 2 = -1 (STI with a surface Dirac cone at the F- 
point) . The lower panel is for the same nanowire pierced by a 
7r-flux tube. To avoid formation of a bound state, the 7r-flux is 
divided into two 7r/2-flux tubes penetrating two neighboring 
plaquettes in the wire center region. 



leading to opening of the gap. Let us explicitly see this. 
The wave function of such a surface electronic state is an 



eigenstate of Eq. (23), which is a plane wave, 



V>0,<£) = e ik * z e ik * R 't', 
The corresponding eigenenergy reads 



(52) 



(53) 



Here, we neglect the anisotropy (B = A). The energy 
spectrum E(k z ) of the surface electronic states is deter- 
mined by imposing (anti-)periodic boundary condition to 
Eq. (52). The spin Berry phase it replaces the periodic 



boundary condition: 

ip(z, 4> + 2ir) = ip(z,(f)), 
by an anti-periodic boundary condition, 



(54) 



(55) 



shifting the allowed values of k^R from integers 
(k<pR = 0, ±1,±2, ••■) to half odd integers [k^R = 
±1/2, ±3/2, • • • ). Importantly, the k^ — and cor- 
respondingly the zero-energy bound state was purged 
from the lowest energy portion of the spectrum by this 
7r-phase shifting (cf. FIG. 3, upper panel). In the 
Aharonov-Bohm geometry, this finite-size energy gap as- 
sociated with the spin Berry phase is compensated by an 
Aharonov-Bohm phase 4>ab , and in some cases we expect 
the energy spectrum closes its gap. 

In the previous section the spin Berry phase 7r for sur- 
face electrons is derived in the presence and absence (due 
to crystal anisotropy) of cylindrical symmetry, that is, 
regardless of fulfilling the spin-to-surface locking condi- 
tion. Here, we confirm numerically the robustness of the 
7T spin Berry phase against cylindrical symmetry break- 
ing, using a tight-binding model for a topological insu- 
lator nanowire. We attempt to show that the amount 
of the "integrated" Berry phase, corresponding to a 27r- 
rotation of the azimuthal angle (j), is precisely 7r. To ver- 
ify this explicitly, we introduce a 7r magnetic flux tube 
piercing the nanowire, as a probe, and investigate the cor- 
responding energy spectrum. Exact cancellation of the 
spin Berry phase and the Aharonov-Bohm phase is con- 
firmed by observing the closing of finite-size energy gap. 

To implement a "shape" in real space such as a 
nanowire we consider in the following a tight-binding ver- 
sion of Eq. (35) on a cubic lattice. We first make the 



following replacement: 



kj — > sin kj , 



(56) 



where j = x,y,z, for the fcj's in Eq. (35), and similarly, 
k) 2(l-cosfcj), (57) 
for fc^'s in M. After this replacement the Hamiltonian 



( 35 1 can be interpreted as a tight-binding Hamiltonian, 



H sd= {(Mo+6M 2 )\x,y,z)(x,y,z\ 

x,y,z 

+ (t x \x + l,y,z)(x,y,z\ + t y \x,y + 1, z){x, y, z\ 
+ t z \x,y,z + l)(x,y,z\ + h.c.)}, (58) 



where 

t, 



.B 

* T~x 0~ x 



M 2 T Z , ty 

M 2 T Z . 



A 



MoT, 



(59) 



Depending on the value of Mq/M 2 , the tight-binding 
Hamiltonian (58) describes either strong/weak topologi- 
cal insulators (STI/WTI) or an ordinary insulator; —4 < 
M /M 2 < and -12 < M /M 2 < -8 STI, -8 < 
M /M 2 < -4 -)• WTI, M /M 2 < -12 and < M /M 2 
— > ordinary insulator. In the following, we will mainly 
focus on the case —4 < M /A1 2 < 0, corresponding to 
STI with a surface Dirac cone at the T-point. So far the 
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V 

FIG. 4. Energy spectrum E(k z ) and a typical shape of the 
(squared amplitude of) the lowest energy wave function in 
the presence of a 7r-flux tube penetrating a series of single 
plaquettes along the axis of nanowire. The cross section is 
rectangular, N x — 12, N y = 16, the anisotropy is B/A — 0.7. 
M0/M2 = — 1. Ah — A = 1. The wave function is depicted 
at k z /n = 0.03. 



Hamiltonian is translationally symmetric. We now re- 
strict the electrons to move only inside the nanowire with 
a rectangular cross section: 1 < x < N x , 1 < y < N y . 
This can be done simply by switching off unnecessary 
hopping amplitudes. After this only k z remains to be 
a good quantum number; we always assume a periodic 
boundary condition in the z-direction. 

We then introduce an Aharonov-Bohm flux <f>AB pierc- 
ing the nanowire. We consider typically the case of 
4>ab — 7T (case of a 7r-flux tube). The simplest way 
to do this is to let a 7r-flux tube penetrate a central 
plaquette of each z-layer, e.g., the plaquette centered 
at (x,y) = (N x /2 + l/2,N y /2 + 1/2) for N x and N y 
being an even integer. The insertion of a flux can be 
achieved by Peierls substitution, e.g., t x —> t x e l ^ AB for 
x = N x /2 N x /2 + 1, y = l,--- ,N v /2. This turns 
out, however, to be not the best solution for us, since 
a 7r-flux tube penetrating a s ingle plaquette involves a 
(zero-energy) bound stateP^H I n the nanowire geometry 
such a bound state appears as the lowest energy helical 
modes bound to the flux tube, propagating in the direc- 
tion opposite to the preformed surface electronic states, 



i.e., electrons localized on the surface of the cylinder. 
Formation of such a pair of counter-propagating modes 
separated only by a finite distance leads naturally to gap 
opening (FIG. 4, upper panel); mixing of the two counter- 
propagating modes (see the wave function, in FIG. 4, 
lower panel) causes level repulsion between the two ini- 
tially gapless states. 

To avoid formation of such a bound state we rather in- 
troduce here a total magnetic flux 7r divided into two 7r/2- 
flux tubes penetrating the two neighboring plaquettes, 
e.g., plaquettes centered at (x,y) — (N x /2 — l/2,N y /2 + 
1/2) and {N x /2 + 1/2, N y /2 + 1/2). With this we could 
see a clear signature of the finite-size gap closing (FIG. 3, 
lower panel). There appears no low-energy bound states 
around a 7r/2-flux tube. In the same figure one can ac- 
tually recognize such bound states pushed up into the 
high-energy spectrum. 

In the above example, by changing the amount of 
Aharonov-Bohm flux tpAB — ^ 27r, 37r, ■ • • , we find an 
even/odd feature (gapped, gapless, gapped, gapless, • • • ) 
in the energy spectrum of surface electronic states. An 
alternative way to verify explicitly that the amount of 
this Berry phase is precisely 7r is to investigate a similar 
even/odd feature due to crystal dislocation lines (results 
not shown here). We have shown previously that the 
electronic states along such a dislocation or equivalently 
a nanowire exhibits a finite size energy gap, manifesting 
the existence of spin Berry phased 

What have we verified in these numerical simulations? 
(Especially in its relation to what we have discussed in 
Sec. II B.) In the example we have presented in this sec- 
tion the cylinder itself is distorted, having a rectangular 
cross section. In addition to such a structural asymmetry, 
we have also taken into account the crystal anisotropy. 
Our data indicate that the global spin Berry phase 7r is 
robust under the coexistence of a structural asymmetry 
and a crystal anisotropy. 



IV. CONCLUSIONS 

We have studied the electronic states on cylindrical 
surfaces of an anisotropic topological insulator. We 
have established an explicit one-to-one correspondence 
between the effective spin in the surface Dirac Hamilto- 
nian and the real spin inherent to the 3D bulk effective 
Hamiltonian. The effective spin on the surface of a topo- 
logical insulator has a property of being constrained on 
its tangential surface, and in particular, when its sur- 
face is warped (e.g., into a cylinder) the effective spin 
feels this change of tangential plane in real space, and 
consequently, the effective spin completes a 2n rotation 
when the reference point travels once around the cylin- 
der. The existence of spin Berry phase 7r is naturally 
understood in this context. However, on the surface of 
a topological insulator with broken cylindrical symmetry 
in the presence of crystal anisotropy, the effective spin 
does not follow locally the tangential plane in real space, 
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i.e., the effective spin can have generally a component 
normal to the surface. This has been shown analytically 
in Sec. II B, using a rather simple model, Whereas, glob- 
ally the spin Berry phase n is robust against anisotropy 
and breaking of the cylindrical symmetry. The latter has 
been verified numerically in Sec. III. 

In the examples we have considered in this paper, the 
"curved" surfaces did not really have a curvature; a cylin- 
drical surface is flat in the proper use of terminology in 
differential geometry (its Riemann curvature is null) . On 
a genuinely curved surface with a finite Riemann curva- 
ture, e.g., on a sphere, the concept of ir spin Berry phase 



may need some modification or a generalization. Naively, 
it is expected to involve a solid angle associated with par- 
allel transport on a sphere. We leave future studies a 
more rigorous discussion on such an issue. 
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